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ABSTRACT 

Based on the modeling of the central emission-line width measured over sub-arcsecond apertures 
with the Hubble Space Telescope, we present stringent upper bounds on the mass of the central 
supermassive black hole, M,, for a sample of 105 nearby galaxies (D < 100 Mpc) spanning a wide 
range of Hubble types (E — Sc) and values of the central stellar velocity dispersion, a c (58 — 419 
km s~ ). For the vast majority of the objects the derived M, upper limits run parallel and above the 
well-known M, — a c relation independently of the galaxy distance, suggesting that our nebular line- 
width measurements trace rather well the nuclear gravitational potential. For values of o~ c between 
90 and 220 km s" 1 the 68% of our upper limits falls immediately above the M, — o~ c relation without 
exceeding the expected M. values by more than a factor 4.1. No systematic trends or offsets are 
observed in this o~ c range as a function of the galaxy Hubble type or with respect to the presence 
of a bar. For 6 of our 12 M, upper limits with o~ c < 90 km s _1 our line- width measurements are 
more sensitive to the stellar contribution to the gravitational potential, either due to the presence of 
a nuclear stellar cluster or because of a greater distance compared to the other galaxies at the low-cr c 
end of the M, — a c relation. Conversely, our M, upper bounds appear to lie closer to the expected 
M. in the most massive elliptical galaxies with values of a c above 220 km s _1 . Such a flattening of 
the M. — a c relation at its high-<r c end would appear consistent with a coevolution of supermassive 
black holes and galaxies driven by dry mergers, although better and more consistent measurements 
for <7 C and K-b&nd luminosity are needed for these kind of objects before systematic effects can be 
ruled out. 

Subject headings: black hole physics, galaxies: kinematics and dynamics, galaxies: structure 



1. INTRODUCTION 

Supermassive black holes (SMBHs) have now been 
discovered in the center of a sufficiently large number 
of nearby galaxies to probe possible links between the 
masses of SMBHs (M.) and the global properties of 
their host galaxies. In fact, i t has emerged that M, 
corre l ates with the luminosity (jKormendv fc Richstone 



19951 iMarconi fc Hunt! I2003h . mass ([Magorrian et ~ 
19981 lHaring fc Rbd 12004ft. s t ellar velocity disper- 



sion (iFerrarese fc Merrittl 120001: iGebhardt et all 12055 



iTremaine et al.ll2002l: iFerrarese fc Fordll2005t) . light con- 
centration dGraham et ah 11200 ID. and gravitational bind- 
ing energy (jAller fc Richstone] uOO?) of the host-galaxy 
spheroidal component, i.e., the entire galaxy in the case 
of elliptical galaxies or the bulge of disk galaxies. In light 
of these findings it is now widely accepted that the mass- 
accretion history of a SMBH is tightly related through 
feedback to the formation and evolution of the host 
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spheroid ("e.g-lSilk fc Re cs 1998; Haehnelt fc Kauffmannl 
[20001: IDi Matteo et'afl feOQS'). some studies having sug- 

;ested a link with the mass of th e dark-matter halo 

Ferraresdl2002T : iPizzella et al.ll2005l ). 

The slope and scatter of all thes e correlations remain 
quite uncertain ([Novak et al.ll2006h . however, in particu- 
lar since they are still based on a limited sample of galax- 
ies with reliable M. that is biased towards early-type sys- 
tems and that is clustered around a rather limited range 
of stellar velocity dispersion (cr) , approximately between 
150 and 250 km s~ . Given the great theoretical interest 
spurred by these findings, there is a pressing need to ac- 
quire better M, statistics, both in terms of the number 
of targets and in terms of broadening the range of parent 
galaxies, in particular towards spiral galaxies. 

Secure M m measurements in external galaxies are tra- 
ditionally obtained through the modeling of the stellar 
and/or gaseous kinematics, most often as derived us- 
ing Hubble Space Telescope (HST) observations in the 
optical domain. The advent of adaptive-optics sys- 
tems working at near-infrared wavelengths has led to 
more st ellar-dynamical measurements of M, fro m the 
ground ([Houghton et alJl2006t HSfowak et al.ll2007t ). Yet, 
such measurements are still quite expensive, not only 
because good-quality measurements of the stellar kine- 
matics in the near-infrared require relatively long ob- 
servations, but also because proper modeling of the 
stellar kinematics in the immediate vicinity of SMBHs 
needs robust constraints on the importance of radial or- 
bits and thus additional large-scale observations, possi- 
bly with integral-field spectroscopy (IValluri et al.l [20041 : 
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ICappellari k McDermidll2005l ). Water -masers have pro- 
vided the most accurate extragalactic M, measurements 
to date, but such gaseous systems a re exceedingly rare 
(jBraatz et al.1 [1991 iGreenhill et all l200l . The mod- 
eling of the nuclear ionized-gas kinemati cs has also 
led to accurate M, measur ements (e.g., IBarth et al.l 
120011: iDalla Bonta et al]|2008t ). and at a relatively cheap 
cost in terms of observation time compared to stellar- 
dynamical M, determ inations (e.g.. iVerolme et "aLll2002t 
IGebhardt et al.l 120031 ). Yet, only an handful of the ob- 
jects targeted by HST turned out to have sufficiently 
regular gas veloc ity fields for the purpose of modeling 
(|Sarzi et al.ll200lD . Thus, unless a large number of galax- 
ies pre-selected to hav e regular nucle ar gas kinematics 
(for instance following IHo et all I2002D is observed with 
HST if and when the Space Telescope Imaging Spectro- 
graph (STIS) is successfully refurbished, it is unlikely 
that the number of galaxies with secure M m measure- 
ments will increase dramatically in the near future. 

The HST Science Archive already contains an un- 
tapped resource that can be used to better constrain the 
black-hole mass budget across the different morphologi- 
cal types of galaxies, which consists of the vast number 
of the STIS spectra from which a central emission-line 
width can be measured. The modeling of this kind of 
data can i ndeed lead to tigh t upper limits on M m , as first 
shown by iSarzi et al.l ()2002D . For this reason we started 
a program aimed at deriving M, upper limits based on 
HST spectra for the largest possible number of galaxies 
and a wide range of morphological types. In this paper 
we present the results based on a sample of 105 nearby 
galaxies for which STIS/G750M spectra in the Ha region 
and measurements of the stellar velocity dispersion were 
available from the HST archive and in the literature, re- 
spectively. Although we will be able only to set an upper 
limit on the M, of our galaxies, the lack of exact mea- 
surements will be compensated for by the large number 
of upper limits when studying SMBH mass-host galaxy 
relationships. 

The paper is organized as follows. In §[2] we describe 
our sample selection and the measurement of central 
emission-line width, before briefly describing our mod- 
eling. We will then present our results and discuss our 
findings in the context of the M, — a relation between 
the SBHM mass and central stellar velocity dispersion of 
the host spheroid in § [3j 

2. DATA COMPILATION AND ANALYSIS 

2.1. Sample selection and data reduction 

In order to assemble the largest possible sample of 
homogeneous measurements of the central emission-line 
width, we queried the HST Science Archive for objects 
with STIS spectra obtained with the G750M grating 
through either the O'.'l x 52" or the 0'.'2 x 52" slit placed 
across the galaxy nucleus, most often at random posi- 
tion angles. This is indeed by far the most common 
configuration in the archive, which always gives us ac- 
cess to the [N n] AA6548, 6583, Ha and [S n] AA6716, 6731 
emission lines. Limiting ourselves to galaxies within 
100 Mpc to minimize the impact of the stellar poten- 
tial on our results, we retrieved data for 177 galaxies 
spanning the whole range of morphological types. When 
available, galactic distances were adopted from mea- 



surements based either on surface-bri ghtness fluctuations 
(iTonrv et al.|[20qot iTomv et al.l l200lD Ceph eid variables 
jFreedman et al.ll2001ft or from lTullvl (|1988l) . Otherwise 
we used the weighted mean recessional velocity corrected 
to the reference frame defined by the microwave back- 
ground radiation from Ide Vaucouleurs et al.l ()1991l . RC3 
hereafter) to derive the distance to our sample galaxies 
by assuming H = 75 kms -1 Mpc^ 1 , £l m — 0.3, and 
ft A = 0.7 . The median distance of the sample galaxies 
is 21.4 Mpc. 

The archival spectra were reduced using IRAJ0 and 
the STIS reduction pipeli ne maintained by t he Space 
Telescope Science Institute (jDressel et al.ll2007l) . The ba- 
sic reduction steps included overscan subtraction, bias 
subtraction, dark subtraction, and flatfield correction. 
Different spectra obtained for the same slit position were 
aligned using IMSHIFT and knowledge of the adopted 
shifts along the slit position. Cosmic ray events and hot 
pixels were remove d using the task LACOS_SPEC by 
Ivan Dokkuml (|2001[ ). Residual bad pixels were corrected 
by means of a linear one-dimensional interpolation us- 
ing the data quality files and stacking individual spectra 
with IMCOMBINE, This allowed to increase the signal- 
to-noise ratio of the resulting spectra.. We performed 
wavelength and flux calibration as well as geometrical 
correction for two-dimensional distortion following the 
standard reduction pipeline and applying the X2D task. 
This task corrected the wavelength scale to the heliocen- 
tric frame too. 

To measure the nuclear emission-line width we gener- 
ally extracted aperture spectra three (0'.'15) and five pix- 
els wide (0"25) centered on the continuum peak, for the 
O'.'l and 0" 2- wide slit cases, respectively. When the spec- 
tra were obtained with a 2-pixel binning read-out mode 
along the spatial direction, we extracted aperture spec- 
tra three pixels wide (0"3) for the O'.'2-wide slit (Table 1). 
The extracted spectra thus consist of the central emis- 
sion convolved with the STIS spatial point-spread func- 
tion (PSF) and sampled over nearly square apertures of 
0"15 x O'.'l, 0'.'25 x 0"2 or 0"3 x 0"2, roughly corresponding 
to a circular aperture with a radius of 0'.'07, 0'.'13, and 
0'.'14, respectively. The wavelength range of our spec- 
tra is either 6482-7054 A or 6295 - 6867 A, depending 
on whether the G750M grating was used at the primary 
or secondary tilt. The instrumental FWHM was 0.87 A 
(omst = 17 km s _1 ) and 1.6 A (cri ns t = 32 km s^ 1 ) for the 
O'.'l and the O'.'2-wide slit, respectively. The atlas of all 
the extracted spectra will be presented in a forthcoming 
paper. 

To place our M m upper limits with the M m — a relation, 
here we consider only galaxies with velocity dispersion 
measurements in the literature, which were available for 
137 objects. We also dropped a further five objects, since 
upon closer inspection they revealed unrelaxed mor- 
phologies. For a number of objects with a sharp central 
surface-brightness profile, the two-dimensional rectifica- 
tion of the spectrum performed during the data reduction 
produced anomalous undulations in t he flux level of con- 
tinuu m of the very central rows (see iKim Quiian o et aLI 
120071 for details) . This introduced also artificial fiuctua- 

5 IRAF is distributed by NOAO, which is operated by AURA 
Inc., under contract with the National Science Foundation 
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tions in the emission-line flux profiles across the nucleus. 
As constraining the concentration of the nebular emission 
is key to our modeling (see ^2.3[) . this problem forced us 
to remove a further eight galaxies from our sample. 



2.2. Measurement of the emission lines 

In order to derive upper l imits on AI m following the 
method of iSarzi et al.l (|2002l . see also £|2.3[) we need to 
measure both the width of the central nebular emission 
and the radial profile of the emission-line flux, so that we 
can gauge both the depth of the potential well and the 
concentration of its gaseous tracer. To side-step the im- 
pact of broad and/or asymmetric emission arising from 
regions much smaller than our resolution limit, we focus 
on the width of the narrow component of the emission 
from forbidden transitions and disregard the broad-line 
emission in our spectra. In the wavelength range of our 
spectra that means measuring the central width and flux 
profile of the [N n] AA6548, 6583 lines since these are usu- 
ally brighter than the [S ii] AA6716, 6731 lines. The [N n] 
doublet also traces the nuclear kinematics better than 
Ha, given that this line could be significantly affected by 
emission from circumnuclear star-forming regi ons (e.g., 
iVerdoes Kleijn et al.ll2(M iCoccato et al.ll2006l ). 

To measure the central width and flux profile of the 
narrow component of the [N ii] lines we fit our spec- 
tra with multiple Gaussians to match both the broad 
and narrow components of all the observed lines, while 
describing the stellar continuum with a low-order poly- 
nomial. A flux ratio of 1:3 was assumed f or the [Nil] 
doub let, as dictated by atomic physics (e.g.. IQsterbrockl 
1989), and in the presence also of [Sn] emission, both 
the [N n] and [S n] doublets were assumed to share a 
common line centroid and width. In most cases only one 
additional Gaussian component was needed in our fits, to 
describe the Ha emission from the broad-line region, al- 
though in many objects also the forbidden [N n] and [S n] 
lines required double-Gaussian profiles. This allowed us 
to describe also lines with Voigt profiles, where tests on 
18 galaxies showed that our narrowest Gaussian compo- 
nent tend to be only slightly broader than the thermal 
component in the Voigt profiles, generally by less than 
20%. The extra Gaussian in the [N n] and [S n] lines was 
also used to isolate the contribution of blue- or redshifted 
wings. To help de blending t he Hq +fNn] region in some 
cases we followed IHo et all (|1997f ) and assigned to both 
the [N n] lines and the narrow Ha emission the line pro- 
file that was predetermined by fitting the [S n] lines. 
The best-fitting Gaussian parameters were derived us- 
ing a non-linear least-squares minimization b ased on the 
robust Levenberg-Mar quardt meth o d (e.g ., iPress et all 
119961) implemented by iMore et all (|1980f >. The actual 
computation was done using the MPFIT algorithm^] im- 
plemented by C. B. Markwardt under the IDL0 environ- 
ment. In objects with conspicuous stellar absorption fea- 
tures that cannot be accounted for by our minimisation 
routine we checked our results against the line- width and 
flux measurements obtained with the GANDALF soft- 

6 The updated version of this code is available on 
|http: / / cow . physics . w ise . edu/ ~craigm/idl/idl . html | 

' Interactive Data Language is distributed by Research System 
Inc. 
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Fig. 1. — Few examples of continuum-subtracted central G750M 
spectra from our spectral atlas illustrating the various fit- 
ting strategies adopted to match Ha, [N Ii] AA6548, 6583 and 
[S n] AA6716, 6731 emission lines. In each panel the red line show 
the overall line blend, whereas the green dashed-dotted lines and 
blue dotted or dashed lines show the adopted narrow and broad 
Gaussian components, respectively. Shown are also the fit residu- 
als, offset for better visibility. For NGC 4548 the nebular emission 
could be match with single Gaussian profiles, for NGC 1961 we 
needed to add a broad Ha component, whereas for NGC 2273 an 
additional broad and blue-shifted component was needed to macth 
the profile of all lines. For NGC 4450 we added two extremely 
broad Gaussians shoulders offset from the center by several thou- 
sands km s -1 , in additio n to a more typical broad Ha component 
(see also IHo ct al. 2000], for a match to the double-peaked profile 
of this LINER 1.9 nucleus). 

war^3 of ISarzi et al.l (|2006l ). adopting either very young 

8 The updated version of this code is available at 
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(300 Myr) or old (10 Gyr) stellar population templates. 
In most cases the measurements agreed within the errors, 
except for IC 342 and NGC 7331 where the Ha absorp- 
tion line is particularly prominent. For these galaxies we 
adopted the GANDALF values. Finally, in defining our 
detection thresholds we compared the amplitude (A) of 
the best-fitting line profile to the noise level (N) in the 
residuals of the continuum fit, adopting as detected only 
those emission lines for which the A/N ratio was larger 
than 3. Figure Q] shows a few sample spectra illustrating 
the various fitting strategies explained above. A more 
detailed description of our emission-line measurements 
for each of our sample galaxy will be presented with our 
spectral atlas. 

In 14 galaxies the nebular emission was too faint for 
it to be detected given the quality of the corresponding 
spectra, and were consequently dropped from our sample. 
Three further galaxies had also to be discarded because 
their line profile could not be well represented as a simple 
sum of Gaussian components. Finally, two galaxies were 
rejected because the radial profile of the flux of the [N n] 
lines was strongly asymmetric and not suitable for mod- 
elling. Table 1 lists the final sample of galaxies analyzed 
in this paper, which comprises of 105 galaxies which 28 
(27%) are classified as ellipticals, 20 (19%) are lenticu- 
lars, and 57 (54%) are spirals. The central velocity dis- 
persion of the ionized-gas component and the size of the 
aperture we measured are also given in Table 1. Prior to 
modeling, the instrumental resolution corresponding to 
the adopted apertures (17 km s^ 1 and 32 km s _1 for the 
O'.'l and 0'.'2 slit widths, respectively) was subtracted in 
quadrature from the observed line- width values to obtain 
the intrinsic gas velocity dispersion. 

Table 2 lists the 74 rejected galaxies. 

2.3. Modeling the central line width 

Assuming that the width of the nuclear emission traces 
the depth of the gravitational well, we can derive strin- 
gent upper bounds on the mass of the SMBHs in our 
sample galaxies thanks to the exquisite spatial resolu- 
tion of HST. Although the stellar contribution to the 
gravitational potential could affect such estimates, the 
fundamental reason for which a lower limit on M. can 
not be set from such a simple measurements is that the 
observed line-broadening may, in principle, be entirely 
due to additional contributions such as non-gravitational 
forces (e.g., gas pressure or magnetic forces). 

In this s t udy w e follow the procedure described in 
iSarzi etaTI (|2002|) . where a detailed description of the 
method can be found. In short, we assume that the ob- 
served line-broadening arises from the motion of ionized- 
gas in a coplanar thin inner disk of unknown inclination, 
where the gas moves in circular and Keplerian orbits 
around the putative SMBHs. For a given radial pro- 
file of the nebular emission, perfectly edge-on disks lead 
to the broadest lines and therefore to a lower estimate 
of M.. Conversely, the M, value needed to explain the 
observed line width diverges to infinity as we approach 
perfectly face-on configurations. Fortunately, such ex- 
treme orientations are statistically rare. Since randomly 
oriented disks have uniformly distributed cosz, it is pos- 
sible to derive la upper and lower limits on M, by sim- 



ply considering models with nearly edge-on (i — 81°, 
cosi = 0.16) and face-on (i = 33°, cosi = 0.84) orien- 
tations, respectively, comprising 68% of the distribution 
of M, values tha t can explain a given line width (e.g., 
ISarzi et al.l[200l . 

In our models we could disregard the effect on the un- 
known position angle of the disk since we extracted our 
spectra in nearly square apertures, and thus assumed 
that the STIS slit was placed along the disk major axis. 

Clearly, for a given disk orientation the concentration 
of the gas tracer impacts heavily on the M, value needed 
to explain a given line width, to the point that no lower 
limit on M, can be set when the gas profile is unresolved. 
This is why the intrinsic emissivity distribution of the 
gaseous disk has to be constrained from the data. As in 
ISarzi et al.l J2002), we assumed an intrinsically Gaussian 
flux profile centered on the stellar nucleus, which makes it 
easier to match the observed flux profile while account- 
ing for instrumental effects. The choice of a Gaussian 
parametrization is also conservative, since cuspier func- 
tions would have led us to estimate smaller M % . For 
instance, adopting an expone ntial profile for the subsam- 
ple of objects studied also bv ISarzi et al.l fj2002f ) leads on 
average to a 10% decrease for the M, estimates. 

In this work we disregarded the contribution of the 
stellar potential, which would lead to tighter upper limits 
on M % - In principle, it is possible to estimate the stel- 
lar mass contribution by deprojecting the stellar surface 
brightness observed in the STIS acquisition images while 
assuming s pherical symmetr y and a constant mass-to- 
light ratio (jSarzi et al.ll2002l ). In practice, however, this 
would only be feasible for a limited number of objects in 
our sample, given the limited quality of the acquisition 
images for most of our sample galaxies, and the pervasive 
presence of dust absorption features, in particular in spi- 
ral host galaxies. Still, the impact of the stellar potential 
is unlikely to change dramatically our M, estimates, in 
particular for the upper limits derived for nearly face- 
on configuration s. For their samp le of nearby galaxies 
(at 8 - 17 Mpc), ISarzi et all (|2002f) found that including 
the stellar mass contribution reduced the median value 
of the M, upper limits by just ~ 12%. For the median 
distance our sample (21.4 Mpc) the stellar mass contri- 
bution to our 33° upper-limits would be ~ 15%. Similar 
considerations would apply to the M, sensitivit y limit 
of our experiment. In the case of the ISarzi et all (|2002f ) 
sample this value was found to be on average 3.9 x 10 6 
Mq, which is well below most of the M. limits derived 
here and comparable to the smallest M. limits obtained 
for the closest objects in our sample. 

To conclude, we note that the range spanned by the 
two M, values delivered by our Keplerian-disk model in- 
cludes also M, that would be estimated under radically 
different assumptions. For instance, the gaseous disk 
model at i — 60° is equivalent to that of an isotropic gas 
sphere in hydrostatic equilibrium (see ISarzi et all 120021 
for details). The M, estimates we derived for i — 33° 
and 81° are listed in Table 1 for the sample galaxies. Al- 
though strictly speaking both values should be regarded 
as upper- limits, we will refer only to the 33° estimates 
as Mm upper limits, hereafter. 



|http: //www.strw .leidenuniv.nl/sauron/ softwar e. html| 
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Fig. 2. — Comparison between our M, upper limits and accurate 
measurements of M, based on the resolved kinematics of gas (filled 
circles), stars (open circles), and water masers (open square) avail- 
able in the literature. Leftward arrows indicate an upper constrain 
rather than a definite value for M. . The upper and lower edges 
of the dotted lines correspond to the M, values that we estimated 
assuming an inclination of i = 33° and 81° for the unresolved 
Keplerian disk, respectively. 

We have determined the la upper and lower confidence 
limits for the M, for randomly orientated disks for 105 
galaxies with measurable spectra and stellar velocity dis- 
persions available in the literature. For 19 galaxies of the 
sample, either M t measurements or upper limits based 
on resolved kinematics were available (Table 1). Fig. [2] 
shows how such measurements compare with our M, lim- 
its, once our values are rescaled accordingly to the dis- 
tances adopted in these studies. Our M. upper-limits 
are consistent within la with such estimates, except for 
NGC 3031 and NGC 4261. Furthermore no systematic 
offset appears when our upper limits are compared with 
similar upper bounds in the literature, rather than defi- 
nite measurements. A particularly complex blend of nar- 
row Ha+ [N n] and broad Ha lines may have biased our 
M, estimates in NGC 3031. 

To place our M, limits on the various versions of the 
M, — a relation, we ap plied to the aperture correction of 
Uorgensen et al.l (|1995| ) to the literature values of stellar 
velocity dispersion in order to obtain the values a c and 
a e and that would have been measured within a circular 
aperture of radius r e /8 and r c , respectively. The effec- 
tive radii r c of the spheroidal components of our sample 
galaxies were taken from various sources in the literature 
(Table 1) except for few disk galaxies for which r c was 
obtain ed from our own photomet ric decomposition (fol- 
lowing [HlndeilAbreueraDIlQQl) of the if -band images 
retrieve d from the archive of the Two Micron All Sky 
Survey (jSkrutskie et al.ll2006l hereafter 2MASS). 

In Figs. [3] and [4] we compare our M, upper limits to 
the M . — a relation, as given by both iFerrarese fc Fordl 
(2005) and lLauer et all ()2007t ). initially to establish the 
validity of our method over a wide range of velocity dis- 
persions. Our upper bounds show a well defined trend 
with both a c and cr c , running closely above the M. — a c 
and M, — a c relations. In the M m — a c plane a Spearman's 



rank coefficient of 0.9 suggests the presence of a correla- 
tion at 9a confidence level whereas a Pearson correlation 
coefficient of 0.8 supports a linear fit to the logaritmic 
data, which returns a slope of 3.43 ±0.21. At first glance 
such a value would imply a shallower trend than found by 
Ferrarese k, Ford ] (|2005l) and a slope closer to that of the 
Lauer et al.l (|2007l) relation, but we need to keep in mind 
that the derived slope could be significantly affected by 
just a few outliers. In particular, for small value of a our 
M m upper-limits could be biased owning to a larger stel- 
lar contribution to the gravitational potential in small 
and distant galaxies. On the other hand, we found that 
our limits appear to parallel particularly well both ver- 
sions of the M t — a relation for 90 < a c , a c < 220 km s~ , 
whereas at lower and higher a a substantial fraction of 
our Mm limits lie either considerably above or almost on 
top of the Mm — a relation, respectively. 

In the following sections we better quantify and inter- 
pret these first considerations. 

3.1. Main trend in the sample 

In the a c interval between 90 and 220 km s _1 our M. 
upper limits appear to correlate particularly well with 
<7 C , paralleling the M, — a c relation. In this a c region 
a value of 0.8 for the Spearman's rank correlation co- 
efficient suggests the presence of a correlation at a 7-a 
confidence level, whereas a Pearson coefficient of 0.8 in- 
dicates that the logarithm values of our M % upper limits 
and <7 C are very likely to be linearly correlated. A linear 
fit in the log a c — log M. plane delivers a best- fitting slope 
of 4.52 ± 0.41 for our M, upper limits, compar ed to the 
4.86 ± 0.43 slope of the IFerrarese fc Fordl (|2005l ) relation, 
with a scatter of 0.39 dex (Figs.[SJi). In the 90 < a c < 220 
km s _1 interval we have 66 M, upper limits, which have 
a median 2.7 times higher than the expected M. value 
(Fig.[5fc). These upper- limits can range from falling short 
of the expected M, values by a factor 3.7 to exceeding 
them by a factor 17.3, although 68% of them actually do 
not top the expected M, values by more than a factor 
4.1 and fall immediately above the M, — a c relation. For 
comparison, by fitting our upper limits in the M. — a c 
plane we obtain a slope of 4.12 ±0.38, very close t o the 
value of 4.13 ± 0.32 found by lLauer et all ff2007h . In 
fact, the parallel trend of our upper limits holds as far 
as a c ~ 300 km s — , with a Spearman's rank coefficient 
of 0.8, a Pearson linear correlation coefficient of 0.8 and 
a linear slope of 3.84 ± 0.28. 

Fig. shows that such a trend holds independent of 
galactic distance - objects as far away as 60 Mpc appear 
to run parallel to the M % — a c relation. In particular, the 
objects at and below 20 Mpc are well distributed. In fact, 
if in this range of a c we perform separate linear regression 
for the three different populations of upper limits with 
D < 30 Mpc, 30 < D < 60 Mpc and 60 < D < 100 Mpc, 
we find slope values that are consistent within the errors, 
namely of 4.05±0.51, 3.51±1.01, 4.52±1.14, respectively. 
This finding shows that the observed nuclear line widths 
do not simply trace an increasingly larger subtended stel- 
lar mass, as galaxies with progressively larger stellar ve- 
locity dispersions are found preferentially at larger dis- 
tances. Instead, the fact that our upper limits scale in 
the same way with a c as precisely-measured M % determi- 
nations indicates that the nuclear emission we measured 




Fig. 3. — Comparison between our M. upper limits and M, — cr c 
relation by Fcrrarcsc & Ford (2005) (thick line) as a function of 
galaxy distance (a) and morphological type (b). The upper and 
lower edges of the dotted lines correspond to M. values estimated 
assuming an inclination of i = 33° and 81° for the unresolved 
Keplerian disk, respectively. Large circles mark galaxies with <r c < 
90 km s _1 that host a nuclear star cluster. The dashed lines show 
to the la (0. 34 d ex) scatt er in M. . Additionally, to follow the 
discussion in A3, H and £13. 2 1 the dot-dashed line shows the 3cr (1.02 
dex) scatter above the M. — a relation whereas the open circles 
point to objects where a nuclear cluster is present. 

arises predominantly in regions of the gravitational po- 
tential that are dominated by the influence of the central 
SMBHs. 

This is not completely unexpected given that a num- 
ber of HST observations revealed that the narrow-line 
regions of active nuclei appear to be quite concentrated 
with scales less than 50 pc, much more s o than th e 
underlying ste llar density profile (see, e.g., iHol 120081 ). 
Most recently, Wal sh et al.l (|2008|) have mapped the be- 
havior of the narrow-line region for galaxies observed 
with multiple-slit STIS observations. They found that 
all galaxies of their sample exhibit a centrally peaked 
surface-brightness profile, with the majority of them 
further showing a marked gradient of the emission-line 
widths within the sphere of influence of the central 
SMBH. The high degree of concentration of the gaseous 
tracer of the gravitational potential is what allows us 
to closely trace the presence of the SMBH even in ob- 
jects where formally its sphere of influence is not re- 
solved. This is similar to the case of the stellar-dynamical 
estimates of M. for M32, which have not significantly 
changed when moving from grou nd- to space-based ob- 
servations (e.g.. IKormen dv 2004, and references therein) 
due to the exceptional concentration of its stellar light 
profile. Actually, that fact that our upper limits run so 
closely to the M, — a c relations also suggests that non- 
gravitational forces do not generally contribute much to 



Fig. 4. — Same as Fig. [3] but now showing the co mparison be- 
tween our M, upper limits and the M, — a c relation of lLauer et al.l 
(pOOTT ) 

the observed line widths (unless for some reason their 
importance scales with <r c ), although their role cannot 
be firmly excluded on a single-case basis. 

Fig. [3]d also shows that in the er c = 90 — 220 km s _1 
range the upper limits derived in galaxies of different 
Hubble types lie neither closer nor further away from the 
M, — a c relation, although only a relatively small num- 
ber of elliptical galaxies falls in this a c interval. Similarly, 
even though only 38% of the spiral and lenticular galax- 
ies in this er c region are unbarred we do not notice any 
systematic tr end with t he pr esence of a bar, unlike what 
was found hv lGrahaml IpOOfih . 

Since our upper limits appear to trace quite closely the 
expected values for M,, we can take advantage of the sig- 
nificant number of galaxies in our sample to understand 
whether the objects that within the present a c range ap- 
pear to show remarkably large or small upper-bounds are 
in fact exceptional. Assuming that our la limits bracket 
symmetrically the expected values of M. and that our 
upper bounds lie 2.7 times above the M, — a c relation, 
with the aid of a Monte Carlo simulation we found that 
16% of our M, upper limits should lie above the M, — a c 
relatio n by more than 3 times its scatter (adopting 0.34 
dex bv lFerrarese fc Ford 2005]), while 8% of them should 
lie below it by more than its scatter. As Figj3] shows, 
only four out of 66 objects in the a c = 90 — 220 km s _1 
range fall that far above the M. — a c relation, with an 
equal number falling below it by more than its scatter. 
Both sets of objects correspond to 6% of the galaxies in 
the considered a c range. 

The previous considerations strongly argue against the 
presence of exceedingly large M, (i.e., above the M, — <r c 
relation by more than 3 times its scatter) in nearby galac- 
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Fig. 5. — Upper panel (a): Comparison between the linear fit to 
the M m upper limits in the range 90 < a c < 220 km s (thick 
line) and the M, — cr c relation by Fcrrarcsc & Ford (2005, thin 
line). Our linear fit to the M, upper limits in the 90 < <r c < 220 
km s — 1 range (shown by the vertical dashed lines) delivers a best- 
fitting slope of 4. 52±0. 41. Lower panels: Distribution of the ratios 
between the measured upper lim its and the values of M , expected 
from the M. — cr c relation by Fcrrarcsc & Ford (2005) for (b) all 
the sample galaxies, (c) the galaxies with 90 < <r c < 220 km s — 1 , 
(d) with cr c < 90 kms" 1 , and (e) with cr c > 220 kms" 1 . The 
median of each distribution is marked by an arrow. 

tic nuclei, and further suggests that galaxies with consid- 
erably smaller M.-budgets (i.e., below the M, — a c rela- 
tion by more than its scatter) may be particularly rare. 
The presence of undermassive SMBHs in field galax- 
ies ha ve been suggested for instance by iVittorini et al.l 
(2005), who argued that in a low galactic-density envi- 
ronment the M. growth may be hampered by the lack of 
gaseous fuel. A p opulation o f und ermassive SMBH was 
discovered also by IVolon tcri (2003) in her simulations of 
the last stages of black-hole mergers, when the binary ex- 
periences a recoil due to asymmetric e mission of gravita - 
tional radiation. According to Fig. 4 of iVolonteril (|2007j ) 
up to 25% of the galaxies with 90 < a c < 220 km s _1 
could contain undermassive SMBHs, with less than 10% 
the expectected M,. Unfortunately, only a handful of 
objects in our sample are massive and close enough (e.g., 
for o~ rmc > 150 km s _1 and D < 20 Mpc) to allow us 
to probe such a low M, regime, where our simulations 
indicate that we should expect only ~ 1% of our upper- 
limits. Furthermore, although the wide range of Hubble 
types and of values for a c spanned by our sample galax- 
ies suggest these are fairly representative of the general 
properties of the nearby population, our sample is al- 
most certainly incomplete, in particular as we probe the 



low-end of the luminosity function where most galaxies 
are found. Our constraints should therefore be regarded 
with caution. 



3.2. The lower end of the M m — er c relation 

At small cr c (< 90 km s _1 ) half of our upper limits sys- 
tematically exceed the expected M. values by 3 times the 
scatter of the M, — a c relation. They are in average larger 
by more than a factor 40 (Fig. [SJi), consiste nt with pre- 
vious works on much more smaller samples (i Sarzi et aLl 
[2001 [SarSi|[200l IVerdoes Kleiin et al.ll2006D . They are 
hosted by NGC 3021, NGC 4245, NGC 5347, NGC 5427, 
NGC 5879, and UGC 1395, which are late-type spirals 
with different degrees of nuclear a ctivity, as measured by 
IHo. Filippenko. fc Sargentl (|1997l . see Table 1). 

We have considered different possibilities related to the 
measurement and modeling of the [N n] A6583 emission 
line to explain the high values of M, found in these 6 
objects. For instance, the presence of broad or asym- 
metric components in our spectra could affect the width 
of the narrow component of the [N n] lines that we mea- 
sured and consequently the M m upper limits giving larger 
masses. A similar bias would be introduced if the ex- 
tent of the flux profile were to be systematically overes- 
timated. Blue asymmetries are observed in the top out- 
liers NGC 5347 and NGC 5427, which are a lso part of 
the sa mple of active galactic nuclei studied bv lRice et al.l 
(2006), who investigated the resolved kinematics of their 
narrow-line region with STIS spectra and first reported 
the presence of blue wings in the [S n] AA6716, 6731 lines. 
In our fits, however, the contribution of such additional 
features was isolated using double Gaussian profiles. As 
regards the flux profile of the [N n] doublet of the small- 
<r c outliers, these are not systematically shallower than 
the other galaxies following the M, — tr c relation in the 
same o~ c range. Therefore, the M. upper limits that we 
have calculated are not biased by either of these effects. 

To explain the largest M m upper limits found at low a c 
values we also considered the impact of the presence of a 
nuclear star cluster (NC), and in general that of a larger 
stellar contribution due to a greater distance. NCs are 
massiv e stellar clusters c oincident with the galaxy photo- 
center (|C6te et al.l(2006l ) that are found in about 75% of 
late- type spiral galaxies (IBoker et al.112002 ) . The ir mean 
effective radius is ~ 3.5 pc ( Boker et al 1 12004]) . small 
enough for them to be complete ly enclosed within the 
central aperture of our spectra. iFerrarese et al.l (2006) 
found a different M, — a c relation for NCs, with similar 
slope but a normalization that is larger by roughly an 
order of magnitude than that the one found for SMBHs. 
The presence of NCs in our low-cr c outliers could there- 
fore explain why they show such high central mass con- 
centrations as indicated by their high M, values. To 
assess the incidence of NCs in the sample galaxies with 
er c < 90 km s _1 we analyzed their surface-brightness ra- 
dial profile obtained with the IRAF task ELLIPSE on the 
STIS acquisition images. For half of the 1ow-ct c outliers 
we could recognize the presence of a NC (NGC 3021, 
NGC 4245, and NGC 5879). On the other hand, we 
could identify a NC only in one (NGC 4212) of the six 
galaxies (17%) which run close to the M, — a c relation 
(IC 342, NGC 2685, NGC 2748, NGC 3982, NGC 4212, 
and NGC 5194). The presence of a NC in the galaxies at 
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the low-cr c end of our sample is shown in Fig. [31 and in 
the case of NGC 3021 and and NGC 5879 it was already 
known (see lScarlata et al.ll2004J: iSeth et al.l l2008. respec- 
tively) . If our limits trace indeed the dynamical signature 
of a NC in thes e nuclei, better d ata and more detailed 
modelling (e.g-. IBarth et al]|2008f) would be required to 
disentagle the contribution of the NC and SMBH to the 
total mass budget. As regards the distance of the low-cr c 
outliers, although we can only rely on distances inferred 
from their recessional velocities it is significant that half 
of them (NGC 5347, NGC 5427, and UGC 1395) are 
found beyond 30 Mpc, whereas all the other 1ow-ct c galax- 
ies are significantly closer, including those for which our 
M. upper limits lie well within 3 times the scatter of 
M. — <7 C relation. 

These findings suggest that part, if not all, of the ex- 
ceedingly large M, values we found at the low-<7 c end of 
the M. — a c relation could be ascribed to a more signif- 
icant stellar contribution to the gravitational potential. 
This is either because of the presence of a nuclear stellar 
cluster (in NGC 3021, NGC 4245, and NGC 5879) or due 
to a larger galactic distance (for NGC 5347, NGC 5427, 
and UGC 1395) than otherwise required to trace the 
Mm — <j c relation at these a c regimes. Therefore, we 
presently do not need to in voke either non-gravitational 
forces (sec Sarzi et al. 2002) or a popu lation of more mas- 
sive SMBHs (see lGreene &: Hq 12006) to explain the ob- 
served flattening of the M m —a c relation at low M, values. 

3.3. The upper end of the M, — <r c relation 

At high <7 C (> 220 km s _1 ) our M. upper limits nicely 
bracket the M. — a c relation (Fig. GJ) and most of them 
are consistent with its scatter (Fig. [3]). In fact, only four 
objects (15%; NGC 2911, NGC 4552, NGC 4594, and 
NGC 5077) fall above the M. — a c relation by more than 
its scatter, with the same number of galaxies falling as 
far below the M. - a c relation (NGC 3998, NGC 4278, 
NGC 6861, and UGC 1841). These outliers do not stand 
out from the rest of the objects with er c > 220 km s _1 
for any obvious property such as morphology, nuclear 
activity or distance. This behavior is suggestive of an 
actual flattening of the high-mass end of M. — <r c rela- 
tion, in particular considering that in the most massive 
and radio-loud galaxies the ionized-gas velocity disper- 
sion can show a si gnificant excess over a pure ly gravita- 
tional model (e.g.. iVerdoes Kleiin et al.ll2006frl . 

The flattening at high-a values is less evident when 
our upper limits are compare d to the shallower M, — a c 
relation of lLauer et al.l ()2007[ ). but is nonetheless present 
upon closer inspection. In particular, excluding objects 
with a c < 90 km s" 1 where the impact of the stellar po- 
tential on our Mm estimates could be more important, 
we found a systematic flattening in the main trend of 
our upper limits as the high-tr e end of the M. — cr plane 
is progressively populated. Specifically, whereas a linear 
fit to objects with a c — 90 — 220 km s _1 yields a slope of 
4.12±0.38 (§3.3), extending this range to 280 km s" 1 , 340 
km s _1 and up to the maximum a c value in our sample 
of 386 km s _1 results in best-fitting values of 3.86 ± 0.29, 

9 In fact, for the two radio-loud galaxies in our sample that were 
also studied by Verdoes Kleijn et al. (NGC 383 and UGC 7115) 
the derived M m upper limits lie above the M. — cr c relation 



3.78 ± 0.27 and 3.56 ± 0.26, respectively. 

This finding would be in agreement with the pre- 
dictions of semi-analytic models for the coevolution of 
SMBHs and galaxies at the highest end of the mass spec- 
trum, whereby gala xies and SMBHs grow m ainly via gas- 
poor, dry mergers (Sch awinski et al.ll2006D . Yet, the be- 
havior of Mm — cr c relation in this regime is still under de- 
bate. In particular, the limited number of galaxies with 
reliable M, measurement in the range M.> 10 9 M Q are 
actually co nsistent with a steepening of the M, — a c re- 
lation (e.g.. lWvithell2lMlDalk Bonta et al.1 |20Q8|). Fur- 
thermore, the cutoff at a c ~ 400 km s" 1 o f the local 
velocity dispersion function ()Seth et al.ll2008ft implies ei- 
ther that SMBHs with M.> 3 x 10 9 M are extremely 
rare or that if they exists their host galaxies should lay 
considerably abo ve the present M, — a c relation. In fact, 
at these regimes lLauer et al.l (|2007T ) argue that the stel- 
lar luminosity L is better suited than er c to trace M # . 
The Mm — cr c relation should steepen at its high-a c end 
if Lauer et al. arguments are correct, since the observed 
<7 C saturates for the most massive of ellipticals while con- 
sidering increasingly large values of L. 

Although our results suggest a flattening of the M, — a c 
relation, we need to keep in mind that systematic effects 
related to the measurement of the bulge properties may 
be significant at the high-cr c end of the M % — cr c plane. In 
particular, the aperture correction for the stellar velocity 
dispersion may be both more important and more uncer- 
tain for the most massive of ellipticals than for smaller 
elliptical and lenticular galaxies. Indeed, giant ellipticals 
tend to have shallower central surface brightness profiles 
than their less massive counterparts, which makes the 
aperture correction more sensitive to the quality and spa- 
tial coverage of the stellar kinematics and to uncertainties 
on the value of the galaxy effective radius, r c . Inciden- 
tally, measurements or r e are also generally less accurate 
for giant ellipticals, due the presence of extended stel- 
lar halos. Ideally, rather than o~ c one would like to have 
a quantity that is more closely connected to the stel- 
lar mass, such as the total if-band luminosity, wh ich is 
also known to relate to M, (|Marconi & Hunt 2003), or a 
direct measurement of a c . Obtaining the Zf-band lumi- 
nosity of our sample galaxies would require much deeper 
images than the available 2MASS data, whereas properly 
measuring er c requires integral-field observations, such 
as those derived in t he case of the SAURON survey 
(|Emsellem et alj|200l . 

3.4. Summary 

Owing to the exquisite spatial resolution of HST and 
to the concentrated character of the ionized-gas emission 
in low-luminosity AGNs, we have been able to set tight 
upper limits on Mm for a sample of 105 nearby galaxies 
(D < 100 Mpc) using STIS/G750M spectra. This sample 
spans a wide range of Hubble types (with 54% of spirals) 
and includes objects with published values for their cen- 
tral stellar velocity dispersion a c . Our main findings are: 

• Independent of the galaxy distance, our M, upper 
limits run parallel and above the M, — er c relation, 
in particular for values of a c between 90 and 220 
km s _1 . The median of the 66 M, upper limits in 
this <7 C regime exceeds the expected M. value by a 
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factor 2.7, with 68% of our upper limits falling im- 
mediately above the M. — a c relation and without 
exceeding the expected M, values by more than a 
factor 4.1. 

• That our nebular line-width measurements trace 
rather well the nuclear gravitational potential, 
makes large samples of M, upper-limit measure- 
ments useful in constraining the frequency of ob- 
jects with exceedingly low or high values of M, 
and in probing the black-hole mass budget across 
the entire Hubble sequence. 

• No systematic trends or offsets are observed in this 
<7 C range as a function of the galaxy Hubble type, 
or with respect to the presence of a bar. Further- 
more, no evidence was found to suggest that the 
largest or smallest M m upper limit in the er c range 
between 90 and 220 kms -1 are actually bracket- 
ing exceptionally high or low values of M m . Thus, 
galaxies with exceedingly high M, budgets must be 
very rare. 

• For cr c values below 90 km s _1 half of our M, up- 
per limits systematically exceed the expected M, 
values by more than a factor 40, consistent with 
previous work on much smaller samples. 

The line-width measurements for such low-cr c out- 
liers are most likely affected by the stellar contribu- 
tion to the gravitational potential, either due to the 
presence of a nuclear stellar cluster or because of a 
greater distance compared to the other galaxies at 
the low-<7 c end of the M m — <j c relation, for which 
our M. upper limits are closer to the expected M, 
values. 

• At the opposite a c end of the M, — er c relation, 
for values of er c above 220 kms -1 , our M. upper 
bounds appear to lie much closer the expected M, 
in the most massive elliptical galaxies, even falling 
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below the M, — er c relation. This flattening is 
less evident when our upper limits a re compared 
with t he shallower M, — er e relation by lLauer et al.1 
(|2007j ). but is nonetheless present upon closer in- 
spection. In particular, excluding objects with 
<7 C < 90 km s~ , we found a systematic flattening in 
the main trend of our upper limits as the high-<7 c 
end of the M, — <7 C plane is progressively populated. 

Although such a flattening of the M, — er c relations 
at its high-cr c end would appear consistent with 
models for the coevolution of supermassive black 
holes and galaxies driven by dry mergers, we 
caution that better and more consistent measure- 
ments for either the X-band luminosity or the 
integrated value of the stellar velocity dispersion 
<7 C within the bulge effective radius r e (both better 
tracers of the bulge mass than er c ) are needed 
before systematic effects can be ruled out. 
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17 


8055 





.25 


X 





2 


66 


X 


52 


91 ± 3 


2.5E8 


9.4E7 






31 


8607 





.30 


X 





2 


15 


X 


10 


79 ± 2 


1.5E7 


6.8E6 






17 







25 


X 


0. 




28 


X 




147 ± 11 


3.1E7 








31 


8607 





.30 


X 





2 


19 


X 


13 


80 ± 4 


3.9E7 


1.0E7 






32 


8236 





.30 







2 


72 


X 


48 


166 ± 10 


4.2E8 


1.0E8 






38 


8236 





.30 







2 


132 


X 


88 


188 ± 7 


6.4E8 


1.2E8 






20 


8228 





.30 







2 


24 


X 


16 


103 ± 6 


4.5E7 


1.2E7 






39 


7361 





25 


X 





2 


21 


X 


16 


49 ± 2 


1.7E7 


5.2E6 






39 


7361 





25 


X 





2 


20 


X 


16 


103 ± 3 


6.7E7 


1.8E7 






40 


7354 





.15 


X 


0. 


1 


10 


X 


7 


426 ± 11 


3.5E8 


5.0E7 


2.7E8(0.3,4.7) 


63 


40 


7403 





25 


X 





2 


25 


X 


20 


205 ± 5 


1.9E8 


3.6E7 






20 


8228 





25 


X 





2 


22 


X 


18 


36 ± 2 


4.4E6 


1.0E6 


< 2.0E7 


64 


20 


8228 





25 


X 


0. 


2 


15 


X 


12 


50 ± 5 


1.3E7 


3.4E6 






41 


7361 





25 


X 





2 


19 


X 


15 


211 ± 3 


2.0E8 


4.0E7 






26 


8607 





.30 


X 





2 


20 


X 


13 


69 ± 9 


2.6E6 


3.4E5 






27 


7361 





25 


X 





2 


18 


X 


15 


139 ± 3 


1.3E8 


3.9E7 


< 5.2E7 


58 


20 


8228 





25 


X 


0. 


2 


4 


X 


3 


59 ± 5 


2.8E6 


4.0E5 






42 


7361 





25 


X 





2 


19 


X 


15 


117 ± 24 


5.0E7 


5.6E6 






31 


8228 





25 


X 





2 


9 


X 


7 


175 ± 3 


4.8E7 


9.1E6 


3.9E7(3.6,4.2) 


65 


38 


8236 





.15 


X 





1 


23 


X 


15 


188 ± 8 


3.8E8 


9.3E7 


5.4E8(4.2,6.6) 


67 


27 


7403 





25 


X 


0. 


2 


19 


X 


16 


333 ± 8 


1.8E8 


5.2E7 






27 


7361 





25 


X 





2 


20 


X 


16 


56 ± 5 


1.7E7 


4.4E6 






27 


7361 





25 




0. 


2 


18 




15 


65 ± 1 


7.4E6 


3.4E6 






43 


8236 





.30 







2 


92 


X 


61 


171 ± 8 


5.5E8 


1.3E8 


< 6.0E8 


43 


27 


7124 





25 




0. 


2 


22 


X 


18 


396 ± 14 


1.4E9 


4.7E8 


1.5E9(0.7, 1.8) 


66 


29 


8607 





.30 







2 


28 




19 


134 ± 5 


1.7E8 


3.4E7 






29 


7361 





25 


X 





2 


18 


X 


14 


98 ± 4 


4.8E7 


1.6E7 


< 7.5E6 


68 


44 


7361 





.2,-, 


X 





2 


37 


X 


29 


138 ± 3 


2.4E8 


6.5E7 







Table 1 - Continued 



Galaxy 


Morph . T . 


Spec. CI. 




Re . 


n^-0 
B 


r e 


an 


Re . 


CT c 


Re . 


rop. 






o 

M m (i = 33 ) 




(i = 81 ) 


M, (low, high) 


Re . 








(Mpc) 




(mag) 


(") 






(km s ) 






(" ) 


( pc ) 


(km s- 1 ) 


(M Q ) 


(M ) 


(M ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


NGC 4459 


SA0+(r) 


T2: 


16.1 


3 


-19.83 


16.7 


V 


10 


179 ± 13 


23 


7361 


. 25 X 0.2 


19 X 16 


216 ± 9 


2.4E8 


3.0E7 


7.0E7(5.7,8.3) 


58 


NGC 4477 


SBO(s):? 


S2 


22.3 


2 


-20.44 


24.4 


V 


10 


154 ± 9 


45 


7361 


n o c: v n o 
u . z o x u . z 


27x22 


110 ± 2 


7.6E7 


1.9E7 






NGC 4486 


E+0-1 pec 


L2 


15.9 


3 


-21.51 


94.9 


B 


2 


349 ± 8 


23 


8666 


0.25 X 0.2 


19 X 15 


630 ± 14 


2.9E9 


8.8E8 


3.8E9(2.8,4.8) 


89 


NGC 4501 


SAb(rs) 


S2 


34.6 


2 


-22.84 


6.0 


V 


10 


152 ± 17 


31 


7361 


0.25 X 0.2 


42 x 33 


99 ± 1 


8.0E7 


1.5E7 






NGC 4507 


(R')SABb(rs) 


S2* 


50.4 


2 


-21.23 


1.9 


K 


5 


156 ± 7 


18 


9143 


0.25 X 0.2 


61 X 49 


54 ± 3 


3.6E7 


7.3E6 






NGC 4526 


SABO(s): 


H 


16.9 


3 


-20.61 


19.7 


V 


10 


208 ± 3 


46 


9163 


0.25 X 0.2 


20 X 16 


320 ± 11 


3.2E8 


6.1E7 






NGC 4548 


SBb(rs) 


L2 


19.2 


3 


-20.63 


35.2 


V 


10 


153 ± 14 


44 


7361 


0.25 X 0.2 


23 X 19 


77 ± 2 


3.8E7 


9.6E6 






NGC 4552 


E0-1 


T2: 


15.4 


3 


-20.36 


29.3 


B 


2 


254 ±11 


29 


8472 


0.25 X 0.2 


19 X 15 


476 ± 18 


1.9E9 


6.9E8 






NGC 4579 


SABb(rs) 


S1.9/L1.9 


24.6 


2 


-21.66 


17.0 


V 


10 


123 ± 16 


31 


7403 


0.25 X 0.2 


30 X 24 


231 ± 2 


2.3E8 


4.3E7 






NGC 4594 


SAa(s) sp 


L2 


9.8 


3 


-21.57 


50.9 


V 


10 


226 ± 6 


38 


7354 


0.15 X 0.1 


7x5 


489 ± 20 


8.6E8 


2.4E8 


1.0E9(0.3,1.7) 


70 


NGC 4596 


SB0+(r) 


L2:: 


29.3 


2 


-20.83 


17.4 


V 


10 


148 ± 14 


47 


7361 


0.25 X 0.2 


35 X 28 


162 ± 11 


2.0E8 


4.0E7 


8. 0E7(4. 0,12.0) 


58 


NGC 4636 


E0-1 


LI. 9 


14.7 


3 


-20.40 


88.5 


B 


2 


178 ± 9 


48 


8472 


0.25 X 0.2 


18 x 14 


254 ± 13 


6.8E8 


2.5E8 






NGC 4698 


SAab(s) 


S2 


17.6 


2 


-19.99 


29.4 


V 


10 


129 ± 9 


28 


7361 


0.25 X 0.2 


21 X 17 


92 ± 2 


8.9E7 


4.3E7 






NGC 4736 


(R)SAab(r) 


L2 


5.2 


3 


-19.83 


45.2 


V 


10 


106 ± 3 


27 


8591 


0.15 X 0.1 


4x2 


90 ± 5 


1.4E7 


3.3E6 






NGC 4800 


SAb(rs) 


H 


13.4 


2 


-18.51 


15.4 


Ft 


11 


111 ± 2 


39 


7361 


0.25 X 0.2 


16 X 13 


95 ± 10 


3.9E7 


3.3E6 






NGC 4826 


(R)SAab(rs) 


T2 


7.5 


3 


-20.55 


5.0 


V 


10 


115 ± 13 


31 


8607 


0.30 X 0.2 


11 X 7 


99 ± 4 


4.2E7 


1.8E7 






NGC 5005 


SABbc(rs) 


LI. 9 


15.7 


2 


-20.79 


5.8 


V 


10 


215 ± 8 


20 


8228 


0.25 X 0.2 


19 X 15 


204 ± 3 


3.4E8 


1.2E8 






NGC 5077 


E3-4 


LI. 9 


41.9 


2 


-20.90 


22.8 


B 


2 


239 ±11 


35 


7354 


0.15 X 0.1 


30 X 20 


397 ± 13 


2.0E9 


4.4E8 


6.8E8(4.0,11.6) 


71 


NGC 5127 


E pec 




67.0 


2 


-21.32 


48.1 


K 


5 


194 ± 5 


49 


8236 






153 ± 7 


5.1E8 


7.7E7 






NGC 5194 


SAbc(s) pec 


S2 


8.5 


2 


-20.99 


47.7 


V 


10 


76 ± 10 


26 


9147 


25 X 2 


10 X 6 & 


36 ± 2 


2.3E6 


4.3E5 






NGC 5248 


SABbc(rs) 


H 


19.2 


2 


-20.78 


0.6 


R 


16 


128 ± 12 


50 


8228 


0.25 X 0.2 


23 X 19 


45 ± 4 


5.1E6 


9.0E5 






NGC 5252 


SO 


SI. 9* 


95.7 


2 


-20.97 


20.7 


J 


12 


183 ± 26 


26 


8055 


0.25 X 0.2 


116 X 93 


197 ± 6 


7.0E8 


1.3E8 






NGC 5283 


SO? 


S2* 


37.0 


2 


-18.72 


8.7 


K 


5 


148 ± 14 


26 


9143 


0.25 X 0.2 


45 X 36 


61 ± 3 


5.9E7 


1.3E7 






NGC 5347 


(R')SBab(ra) 


S2* 


34.6 


2 


— 19.59 


20.5 


V 


10 


70 ± 13 


26 


9143 


0.25 X 0.2 


42 X 33 


53 ± 3 


4.6E7 


6.5E6 






NGC 5427 


SAc(s)pec 


S2* 


38.9 


2 


-21.22 


29.7 


V 


10 


70 ± 12 


26 


9143 


0.25 X 0.2 


47 X 38 


78 ± 1 


8.1E7 


2.0E7 






NGC 5490 


E 




70.0 


2 


-21.34 


18.5 


B 


2 


292 ± 28 


29 


8236 


0.30 X 0.2 


102 X 68 


278 ± 21 


1.3E9 


2.6E8 






NGC 5643 


SABc(rs) 


S2* 


18.6 


2 


-21.11 


2.6 


H 


8 


97 ± 1 


51 


9143 


0.25 X 0.2 


22 X 18 


95 ± 3 


4.7E7 


3.1E6 






NGC 5695 


S? 


S2* 


58.5 


2 


-20.46 


7.1 


K 


5 


148 ± 2 


26 


9143 


0.25 X 0.2 


71 X 57 


67 ± 2 


2.2E8 


5.1E7 






NGC 5728 


SABa(r): 


S2* 


40.3 


2 


-21.37 


56.0 


V 


10 


210 ± 44 


52 


8123 


0.15 X 0.1 


29 X 19 


121 ± 7 


2.4E8 


6.2E7 






NGC 5879 


SAbc(rs):? 


T2/L2 


11.4 


2 


-18.88 


3.0 


B 


13 


59 ± 9 


13 


8607 


0.30 X 0.2 


17 X 11 


73 ± 4 


8.5E6 


2.4E6 






NGC 6300 


SBb(rs) 


S2* 


15.2 


2 


-20.71 


6.0 


K 


6 


94 ± 5 


18 


9143 


0.25 X 0.2 


18 X 15 


75 ± 5 


2.4E7 


9.4E6 






NGC 6500 


SAab: 


L2 


39.1 


2 


-20.50 


2.1 


K 


5 


230 ± 6 


27 


7354 


0.15 X 0.1 


28 X 19 


168 ± 19 


4.3E8 


1.3E8 






NGC 6861 


SA0-(s): 




28.1 


3 


-20.32 


14.9 


K 


5 


419 ± 15 


53 


9163 


0.30 X 0.2 


41 X 27 


623 ± 61 


1.5E9 


3.6E8 






NGC 6951 


SABbc(rs) 


S2 


17.1 


2 


-20.45 


5.4 


V 


10 


104 ± 11 


24 


8228 


0.25 X 0.2 


21 X 17 


72 ± 1 


1.4E7 


5.9E6 






NGC 7052 


E 




61.7 


2 


-21.26 


34.7 


J 


12 


239 ± 23 


54 


8236 


0.15 X 0.1 


45 X 30 


235 ± 14 


1.2E9 


3.1E8 


3.7E8(2.2,6.3) 


72 


NGC 7331 


SAb(s) 


T2 


13.1 


3 


-21.21 


48.6 


V 


10 


126 ± 4 


34 


8228 


0.25 X 0.2 


16 X 13 


204 ± 26 


1.7E8 


7.4E7 






NGC 7626 


E pec: 


L2:: 


40.8 


2 


-20.99 


38.7 


B 


2 


253 ±11 


38 


8236 


0.30 X 0.2 


59 X 40 


279 ± 16 


8.7E8 


1.9E8 






NGC 7682 


SBab(r) 


S2* 


63.5 


2 


-20.34 


9.4 


K 


5 


122 ± 17 


26 


9143 


0.25 X 0.2 


77 X 62 


112 ± 3 


4.8E8 


1.9E8 






UGC 1214 


(R)SAB0+(rs): 


S2* 


64.2 


2 


-20.47 


37.0 


I 


9 


116 ± 15 


26 


9143 


0.25 X 0.2 


78 X 62 


76 ± 1 


8.9E7 


3.9E7 






UGC 1395 


SAb(rs) 


SI. 9* 


65.2 


2 


-20.21 


2.4 


K 


5 


70 ± 6 


18 


9143 


0.25 X 0.2 


79 X 63 


24 ± 3 


1.4E7 


3.8E6 






UGC 1841 


E 




80.3 


2 


-21.28 








348 ±29 


55 


8236 


0.30 X 0.2 


117 X 78 


330 ± 26 


4.9E8 


2.0E8 






UGC 7115 


E 




94.5 


2 


-20.70 


10.9 


B 


2 


205 ± 38 


55 


8236 


0.30 X 0.2 


137 X 92 


318 ± 14 


3.0E9 


3.9E8 






UGC 12064 


SO - : 




64.4 


2 


-20.19 


2.7 


K 


5 


281 ± 19 


55 


8236 


0.15 X 0.1 


47 X 31 


426 ± 21 


2.5E9 


2.3E8 







a 

CD 



c 
p 

o 

B 

r s. 
r s. 
a 
r s< 

O 



O 

td 

o 
w 

K 

o 

ST 
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Notes . — Col. ( 1) : G a laxy name. Col. (2) : Morphological type from RC3. Col. (3) : Nuclear spectral class from 
IHo. Filippenko. & Sargenti fl997h . where H = H II nucleus, L = LINER, S = Scyfert , T = transiton object (LINER/HII), 1 
— type 1, 2 — type 2, and a fractional number between 1 and 2 denotes various intermediate types; uncertain and highly uncertain 
classifications arc followed by a single and double colon, respectively. The nuclear spectral class of galaxies marked with * is from 
NASA/IPAC Extragalactic Database (NED). Col. (4): Distance. Col. (5): Reference for col. 4. All the distances were taken from 
literature (see attached list), except those we obtained from V3K, the weighted mean recessional velocity corrected to the reference 
frame of the microwave background radiation given in RC3. These were derived as Vsk/Hq with Ho — 75 km s" 1 Mpc — . Col. (6): 
Absolute corrected B magnitude derived from (RC3) with the adopted distance. Col. (7): Effective radius of the spheroidal 
component. 

Col. (8): Band in wich the effective radius were measured. Col. (9): Reference for col. (7). All the effective radii were taken from 
literature (see attache d list), except for those we measured by a photometric decomposition of the K-b&nd images available in the 
2MASS science archive Skrutskic ct al. (2006). Col. (10): Central velocity dispersion of the stellar component within r e /8. Col. (11): 
Reference for the measured stellar v elocity dispersion and c orresponding size of the central aperture from which we calculated the 
value given in col. (10) by following J0rgcnscn et alj (I1995T ■ We did not apply any aperture correction to the measured stellar 
velocity dispersions of NGC 2748, NGC 3982, and UGC 1841, because no information about the size of the aperture was available. 
Col. (12): HST proposal number under which was obtained the STIS/G750M spectrum from which we measured the central velocity 
dispersion of the ionized gas. Col. (13): Size of the central aperture where we measured the velocity dispersion of the ionized gas. 
Col. (14): Physical size of the central aperture where we measured the velocity dispersion of the ionized gas. Col. (15): Central 
velocity dispersion of the ionized-gas component within the aperture in Col. (13). This is the intrinsic velocity dispersion obtained 
from the observed one by subtracting the instrumental velocity dispersion. Col. (16): Mm upper limit for a Kcplcrian disk model 
assuming i — 33°. Col. (17): Af, upper limit for i — 81°. Col. (18): Mass (and confidence interval) of the SMBH derived from 
modeling based on the resolved kinematics. The M m of NGC 3227 and NGC 4258 were obtained by studying the dynamics of stars 
and water masers, res pectiv e ly. Th e ioni zed-gas dynamics was u sed fo r all t he remaining galaxi es. Co l. (19): Reference for col (18). 
References. — Qfc jTullvl 11988D: tl) Ide Vaucouleurs et alJ ( | 199ll); p|) ITonrv et alj I2001T): (I) fFreedman et aU l|20"oiT): 0) 
Skrutskic et al. (2006); (6) Kassin ct al. (2006); (7) do Souza ct al. (2004); (8) Laurikaincn ct al. (2004); (9) Xanthopoulos 
11996 1 ) ; 1 101) IBaggett et alj 119981); dill) lAndredakis fc Sanders! 119941 1 : U2 I ) | Marconl &; Hunti 120031); 1131) lFalc6n-Barroso et ahl 
120021); 1141) ILaurikainen et all 12005 | ); 1151) IFisher et all 119951): 119) [Scarlata et al l 1200411 ; 11171) ITerlevich et alj U99oT) ; 
1 181) IGarcia^Rissrnann et alj 120051); 119) [wegner et alj I2003B: feet) IBatcheldor et alj (l2005l);12ll) IDavies et all (119871); (1221) 
Sin.:. 1 1 & Prugniel (1997a); (23) Bcrnardi ct al. (2002); (24) Schechtcr (1983); (25) Bonfanti et al. (1995); (26) Nelson & Whittle 

(1995) ; (27) Barth ct al. (2002); (28) Corsini ct al. (1999); (29) Simien & Prugniel (1997b); (30) Simien & Prugniel (1997c); 
hi!) IrTc raudcau fc 'Simienl ( 1199a): 132n Idi Nella et al l 11995): (33) Hcra udeau et alj 11999)) : 1341) IVega Beltran et ali 120011); 1351) 

ICarollo et all 119931); 139) ISimien fc Prugniell 119981); 1371) IWhitmore et alj 119791); 1381) ISmith et I all 120001); 1391) ISarzi et alj 
120021); 1401) (FisheJ 119971); J41I) ISimien &: Prugniel! 120021) ; 14 2) IFalcon-Barroso et all 120061); 1431) IVerdoes Kleiin et all ]2002l) : 
(44) Whitmorc & Kirshncr (1981); (45) Jarvis et al. (1988); (46) Proctor & Sansom (2002); (47) Bcttoni & Gallctta (1997); (48) 
Bender ct al. (1994); (49) Nocl-Storr ct al. (2007); (50) Dumas ct al. (2007); (51) Gu ct al. (2006); (52) Wagner & Appcnzcller 
119881) ; 1541) [van den Bosch fc van der Marell 1199a): 1551) IBalcells et all 119951); 1551) ITonrv fc DavS 119811); ^571) [Atkinson et all 
I2005I): 1581) ISarzi et all (120011): 1591) IDevereux et all 120031): l60f) IDavies et al l 120061): 16lT) IBarth et all 12001J): ieT) 
Pastorini et al. (2007); (63) dc Francesco et al. (2006); (64) Marconi et al. (2003); (65) Mivoshi et al. (1995); (67) Fcrrarese et al. 

(1996) ; (68) Coccato et al. (2006); (69) Macchetto ct al. (1997); (70) Kormcndy (1988); (71) dc Francesco ct al. (2008); (72) 
van der Marcl & van den Bosc h (1998); (66) Bower ct al. (1998); (53) Koprolin & Zcilinger (2000). 
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TABLE 2 

Properties of the galaxies rejected from the main sample 



Galaxy 


Morp. T. 


Prop. 


Rej. 




(RC3) 






IC 5096 


Sbc sp 


9046 


ns 


NGC 0134 


SABbc(s) 


8228 


ns 


NGC 0157 


SABbc(rs) 


8228 


ns 


NGC 0255 


SABbc(rs) 


8228 


ns 


NGC 0278 


SABb(rs) 


7361 


s 


NGC 0520 


pec 


8669 


i 


NGC 1097 


SBb(s) 


9782 


ds 


NGC 1255 


SABbc(rs) 


8228 


n 


NGC 1300 


SBbc(rs) 


8228 


P 


NGC 1832 


SBbc(r) 


8228 


ns 


NGC 2623 


pec 


8669 


i 


NGC 2654 


SBab sp: 


9046 


s 


NGC 2892 


E+ pec: 


8236 


s 


NGC 2976 


SAc pec 


8591 


s 


NGC 3003 


Sbc? 


8228 


s 


NGC 3049 


SBab(rs) 


7513 


s 


NGC 3067 


SABab(s) 


8596 


f 


NGC 3162 


SABbc(rs) 


8228 


r 


NGC 3254 


SAbc(s) 


8228 


n 


NGC 3256 


pec 


8669 


i 


NGC 3259 


SABbc(rs) 


8228 


s 


NGC 3379 


El 


8589 


f 


NGC 3403 


SAbc: 


8228 


ns 


NGC 3489 


SAB0+(rs) 


7361 


P 


NGC 3516 


(R)SBO(s) 


8055 


d 


NGC 3521 


SABbc(rs) 


8228 


ns 


NGC 3684 
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Notes. — Col.(l): Galaxy name. Col. (2): Morphological type from RC3. Col. (3): HST proposal number under which was obtained 
the STIS/G750M spectrum. Col. (4): Reason of rejection, where d — problem in dcblcnding the emission lines, i — interacting galaxy, 
f — irregular or strongly asymmetric radial profile of the flux of the [N II] emission linca n — faint or absent emission lines, p — 
double-peaked emission lines, r — unsuccessful two-dimensional rectification of the spectrum, s — no stellar velocity dispersion 
available in literature. 



